Abstract. We have analyzed the available polarization 
surveys of the Galactic emission to estimate to what ex- 
tent it may be a serious hindrance to forthcoming ex- 
periments aimed at detecting the polarized component of 
Cosmic Microwave Background (CMB) anisotropics. Re- 
gions were identified for which independent data consis- 
tently indicate that Faraday depolarization may be small. 
The power spectrum of the polarized emission, in terms 
of antenna temperature, was found to be described by 
■ ^ ; Cf ~ (1.2 ± 0.8) • 10-9 • (€/450)-i-«±°-3 • (r//2.4GHz)-5-8 
(^ . K^, from arcminute to degree scales. Data on larger an- 
^^ I gular scales {i < 100) indicate a steeper slope ~ £~^. We 
Cn , conclude that polarized Galactic emission is unlikely to be 
5h ' a serious limitation to CMB polarization measurements at 
^ ,: the highest frequencies of the MAP and Planck-LFI in- 
^5 ' struments, at least for ^ > 50 and standard cosmological 
fy-\ . models. The weak correlation between polarization and 
p\J ' total power and the low polarization degree of radio emis- 
, sion close to the Galactic plane is interpreted as due to 
^N ' large contributions to the observed intensity from unpo- 
.^ , larized sources, primarily strong Hil regions, concentrated 
*y-s ' on the Galactic plane. Thus estimates of the power spec- 
,—1 , trum of total intensity at low Galactic latitudes are not 
0^ ' representative of the spatial distribution of Galactic emis- 
^^ ,: sion far from the plane. Both total power and polarized 
3^ ' emissions show highly significant deviations from a Gaus- 
■»..., sian distribution. 

'J^. Key words: Polarization - ISM: structure - Galaxy: gen- 
O ' eral - cosmology: cosmic microwave background - radio 
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1. Introduction 

Several ongoing or planned experiments (see Staggs 
et al.l999 for a recent review) are designed to reach the 
sensitivities required to measure the expected linear po- 
larization of the Cosmic Microwave Background (CMB). 
The forthcoming space missions Planck and MAP, 
aimed at obtaining full sky high sensitivity and high res- 
olution maps (FWHM of about 56', 41', 28', 21', and 13' 
for MAP channels at 22, 30, 40, 60, and 90 GHz, re- 
spectively; of 33', 23', 14', 10' for Planck "radiomet- 
ric" channels at 30, 44, 70, and 100 GHz; of 10'. 7, 8', 
5'. 5, 5', 5', 5' for Planck "bolometric" channels at 100, 
143, 217, 353, 545, and 857 GHz, respectively) of CMB 
anisotropics will also probe the CMB polarization fluctu- 
ations (Mandolesi_e^alJ^98^^uge^ al. 1998; MAP web- 
page: |http: //map. gsfc.nasa.gov/; Planck webpage: 



tittp: //astro . estec . esa.nl/SA-general/Pro"jects/ - 
Planck/). We recall that the CMB radiation brightness 
peaks at about 160 GHz. 

The current design of instruments for the Planck 
mission (the third Medium-size mission of ESA's Horizon 
2000 Scientific Programme) provides good sensitivity to 
polarization at all LFI (Low Frequency Instrument) fre- 
quencies (30-100 GHz) as well as at three HFI (High Fre- 
quency Instrument) frequencies (143, 217 and 545 GHz). 
The NASA's MIDEX class mission MAP has also polar- 
ization sensitivity in all channels. 

While there is a very strong scientific case for CMB 
polarization measurements (cf., e.g., Zaldarriaga 1998 and 
references therein) , they are very challenging both because 
of the weakness of the signal and because of the contam- 
ination by foregrounds that may be more polarized than 
the CMB. 

Our knowledge of polarized foreground components is 
very meager (see Davies & Wilkinson 1999 for a recent 
review). In this paper we present a preliminary investi- 
gation of the power spectrum of the polarized Galactic 
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synchrotron emission, the likely dominant foreground con- 
tribution at microwave frequencies where CMB is domi- 
nating, at least at intermediate to large angular scales 
(Tegmark et al. 2000). When this work was approaching 
completion we learned that a similar analysis was carried 
out by Tucci et al. (2000). We improve on their results 
by taking into account, in addition to the Parkes survey 
(Duncan et al. 1995, 1997; hereafter D97) discussed by 
them, the more recent Effelsberg surveys at 2.7 GHz (Dun- 
can et al. 1999; D99) and at 1.4 GHz, covering areas up 
to ±20° of Galactic latitude (Uyaniker et al. 1998, 1999; 
U99), as well as the Leiden surveys (Brouw & Spoelstra 
1976; BS76). We also discuss in some detail the effect of 
the emission from Hii regions in the Galactic plane and of 
the Faraday depolarization. 

The Galactic polarized thermal dust emission has been 
modelled by Prunet et al. (1998); based on their results, 
we may expect that fluctuations of polarized dust emission 
prevail over those of polarized synchrotron emission above 
100-150 GHz. De Zotti et al. (1999) discussed polariza- 
tion fluctuations due to extragalactic sources. Additional 
polarized contributions are expected from magneto-dipole 
emission or rotational emission of dust grains (Draine & 
Lazarian 1999; Lazarian & Draine 2000) and scattered 
free- free emission (Keating et al. 1998). A multifrequency 
Wiener flltering method to detect CMB polarization in the 
presence of polarized foregrounds has been worked out by 
Bouchet et al. (1999). 

2. Polarization surveys 

Linear polarization observations extending up to high 
Galactic latitudes and carried out at several frequencies, 
from 408 to 1411 MHz, have been presented by BS76 (Lei- 
den surveys). The half power beamwidths varied with in- 
creasing frequency from 2.3° to 0.6°. Unfortunately these 
and the other surveys discussed by Spoelstra (1984) are 
undersampled so that proper smoothing to the largest 
beamwidth, as required to combine data at different fre- 
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quencies, is not possible. Thus, estimates of differential 
polarization and differential Faraday rotation across the 
beam cannot be made. We have limited our analysis of 
these data to patches of the sky with better than average 
sampling. 

High resolution polarimetric surveys of strips around 
the Galactic plane at 2.4 and 2.695 GHz, respectively, 
have recently been pubhshed by D97 and D99. U99 car- 
ried out a surface brightness and polarization survey of 
four fields at medium Galactic latitudes (up to 6 = 20°), 
at 1.4 GHz. These data, together with details about the 
instrumental capabilities, can be found on the WEB site 



tittp : //www . mpif r-bonn . mpg . de/survey . html 



D97 covered 127° of Galactic longitude (238° < 1 < 5°) 
out to at least b = ±5° (for 340° < 1 < 5° the survey 
extends to 5 = +7°, for 340° < 1 < 352° to 6 = -7°, 
and for 240° < 1 < 270° to b = -8°) with an an- 
gular resolution of 10.4'. The surveyed area amounts to 
~ 1413 square degrees. The nominal rms noise in total 
power is 17 mJy/beam area (8 mK), and 11 mJy/beam 
area (5.3 mK) for polarization; however a lower rms noise 
of 11 mJy/beam area (5.3 mK) for total intensity and 

6 mJy/beam area (2.9 mK) for polarization has been 
achieved over ~ 43% of the total area. The center fre- 
quency is 2.417 GHz and the bandwidth about 145 MHz. 
Values of the Stokes Q and U parameters and of the total 
power are given every 4' on a rectangular grid, in units of 
mJy/beam area; the conversion factor to brightness tem- 
perature is 1 mJy/beam area = 0.48 mK (Duncan, private 
communication) . 

Polarimetric data from the Effelsberg 2.695 GHz sur- 
vey with half-power beamwidth of 4.3' in the first Galactic 
quadrant were reported by D99. Maps at a resolution of 
5.1' of the Stokes Q and U emission components, covering 
the region 74° > 1 > 4°. 9, |6| < 5°, with a rms noise of 
2.5 mJy/beam area (corresponding to 9 mK) were con- 
structed. 

Of particular interest for our purposes is the contin- 
uum and polarization survey at 1.4 GHz, also carried out 
with the Effelsberg 100-m telescope, at medium Galac- 
tic latitudes (up to 6 = ±20°). Four areas were observed 
(with the one in the Cygnus region split in two parts), 
totaling about 1050 sq. deg.: one in the first Galactic 
quadrant (45° > 1 > 55°, 4° > 6 > 20°); the north- 
ern (65° > 1 > 95°, 5° > 6 > 15°) and the southern 
(70° > 1 > 100°, -15° > 5 > -5°) parts of the Cygnus re- 
gion; the highly polarized "fan region" (140° > 1 > 153°, 
3.5° > 6 > 10°); the anticentre region (190° > 1 > 210°, 
3.8° > b > 15°). The rms noise is about 15 mK (about 

7 mJy/beam area) for total intensity and about 8 mK in 
linear polarization; the angular resolution is 9.35'. 

From the MPIfR survey Web-site mentioned above it 
is possible to download both "background" and "source" 
total power maps, which correspond to large and small- 
scale components of the total intensity, as well as polariza- 
tion maps. To derive the total Galactic emission we have 



0= I 0.15 
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Fig. 1. Total versus polarized intensity in the 300° < 1 < 
320° region from the D97 data at 2.4 GHz, showing no 
correlation at all. 

subtracted from the sum of "background" and "source" 
maps the isotropic component due to the CMB and to the 
contribution of unresolved extragalactic sources. We have 
adopted, for this component, Ta, isot = 2.9 K. 

3. Analysis of survey data 

3.1. Depolarization 

A comparison of total power and polarized emission maps 
(Duncan et al. 1995, 1997; Uyaniker et al. 1998, 1999) 
shows little correlation. While the total intensity clearly 
peaks on the Galactic Plane (apart from a number of spurs 
and loops extending to high Galactic latitude), the polar- 
ized intensity is much more uniformly distributed. Also, 
many sources which are very intense in total power are not 
seen in polarized emission and, conversely, bright regions 
of extended polarization do not appear to be connected 
with sources of total-power emission (D97, U99). This is 
seen in Fig. 1, showing, as an example, the polarization in- 
tensity \/Q^ + V^ versus the total intensity for the Parkes 
survey data in the region 300° < 1 < 320° and |&| < 5°. 

Two main factors may contribute to this situation. On 
one side, many bright structures on the Galactic plane 
are (unpolarized) Hii regions and significant thermal ra- 
dio emission is also expected between and above bright 
Hii regions in our Galaxy: Duncan et al. (1995) estimate 
a thermal flux level of the order of 100 mJy per beam area 
at 2.4 GHz, comparable to the level of the radio contin- 
uum often observed near these regions. The large thermal 
contributions to the Galactic emission, which are concen- 
trated close to the Galactic plane, obviously make it very 
unlikely that the total intensity power spectra for these 
regions can be representative of the synchrotron power 
spectra at high Galactic latitudes. 

On the other side, differential Faraday rotation or vari- 
ations of the magnetic field orientation may strongly de- 
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polarize the emission from distant regions of the Galaxy 
(Burn 1966, Gardner & Whiteoak 1966; for a recent, 
detailed discussion of depolarization mechanisms, see 
Sokoloff et al. 1998) so that only the polarized emission 
of relatively local origin can be observed. The two factors 
may act together: variations in the density of thermally 
emitting electrons may lead to a large enough differential 
Faraday rotation to produce substantial depolarization; in 
addition, the magnetic field may be tangled by turbulent 
motions of the ionized gas, leading to further depolariza- 
tion. 

The Faraday rotation of the polarization position angle 
of a linearly polarized wave at a wavelength A traversing 
an ionized medium with electron density Ug and a regular 
magnetic field B is given by: 
0(A) = RMA^ rad, (1) 

where the rotation measure RM is the linc-of-sight in- 
tegral 



KM = 



27rTO2c^ 



=B • dl 



== 810 / ne{cm-^)Bn{nG) ■ dl{kpc) radm"2_ (2) 

Since (j> scales as j^~^, the Faraday rotation is likely 
irrelevant at microwaves, where MAP and Planck in- 
struments will operate, while it may strongly distort the 
power spectrum of polarized emission at the decimetric 
wavelengths considered here, on one side by wiping out 
polarized emission from far regions (due to differential 
rotation within the beamwidth or the bandwidth of ob- 
servations) and, on the other side, introducing structure 
on a variety of scales (larger than the beamwidth) due 
to fluctuations in the thermal electron density distribu- 
tion and/or magnetic field variations (in strength and/or 
direction) along the line of sight. 

If the synchrotron emission arises throughout the 
depth of the Faraday rotating medium, the polarization 
degree is reduced from the intrinsic value Pq to (Burn 
1966): 



P{X) = Po 



sin( 



(3) 



For (p — 1 rad, corresponding to RM(radm ^) = 21.8, 
64.1, and 80.8 for v ^ 1.4, 2.4, and 2.7 GHz, respec- 
tively, the depolarization amounts to about 16%. Spoel- 
stra (1984), from his analysis of multifrequency data on 
linear polarization, found that the distribution of RM over 
the sky shows a complex pattern with typical values (in 
the region covered by Leiden surveys) of 8radm^^. As 
pointed out by the referee, however, Spoelstra's estimates 
of RMs are actually lower limits because there may be 
field reversals along the line of sight, his surveys have dif- 
ferent resolutions and differential Faraday depolarization 
can lead to small observed RMs (see Sokoloff et al. 1998). 
Accurate estimates of the amount of depolarization were 
obtained, for the Leiden survey at 820 MHz, by Berkhui- 
jsen (1971) who determined the percentage polarization 



Table 1. Polarization degrees of non-thermal emission at 
1.4 GHz for several U99 fields (Bcrkhuijsen, private com- 
munication) . 



1 (dcg) 


b (deg) 


P(%) 


151.5 


+6 


20 


146 


+7 


4 


50 


+10 


3 


78 


+10 


19 


86 


+10 


12 


80 


-10 


9 


92 


-10 


10 



using the absolutely calibrated total power survey at that 
frequency. The polarization percentages given in that work 
should be scaled by 0.8 (Bcrkhuijsen 1975) and are not 
corrected for the thermal contribution from ionized gas, 
which decreases the polarization degree. 

The implied RMs are relatively small over the surveyed 
area, particularly at high Galactic latitudes. 

The Galactic RM distribution is also probed by Fara- 
day rotation measurements towards pulsars and extra- 
galactic radio sources. Since pulsar distances can often 
be independently derived, pulsar data also provide infor- 
mation on the variation of RM along the line of sight; 
also, they appear to have no intrinsic Faraday rotation 
and hence their observed RM arises entirely along the path 
to the observer. Extragalactic sources can provide infor- 
mation on the Galactic medium out to large distances, 
beyond those where pulsars are found; on the other hand, 
extragalactic sources may have their own Faraday rotation 
which adds to the Galactic contribution. 

A catalogue of known pulsars, including values of RM, 
has been published by Taylor et al. (1993); an updated 
version is available on the WEB (see Appendix of Tay- 
lor et al. 1993). Additional RMs have been published by 
Manchester & Johnston (1995), Navarro et al. (1997) and 
Han et al. (1999); on the whole, we have collected RMs 
for 318 pulsars. 

Rotation measures for 674 extragalactic sources have 
been catalogued by Brotcn et al. (1988). Additional RMs 
have been published by Clegg et al. (1992). 

An analysis of the RM of pulsars and extragalactic 
sources located in the sky areas covered by polarization 
surveys considered here reveals a number of regions where 
RMs can produce only a small Faraday depolarization 
(</) < Irad). 

One of these is the area 140° < 1 < 153°, 3°. 7 <b< 
10°, surveyed by U99. This partly covers the highly polar- 
ized region referred to as the "fan region" where the rota- 
tion measures have long been known to be small (Bingham 
& Shakeshaft 1967) and the magnetic field direction has 
to be basically perpendicular to the line of sight. 

Bcrkhuijsen (private communication) has estimated 
the observed polarization degree of non-thermal emission 
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Fig. 2. Distributions of polarized intensities for the U99 
regions centered at 1 = 50° (solid line), 80° (6 > 8°; short 
dashes), 146°. 4 (dotted line), and 200° (long dashes). 

for several U99 fields discussed in this paper (see Table |l|). 
The total synchrotron emission was derived subtracting 
from the observed intensities the contributions from ex- 
tragalactic radio sources and the CMB (3 K at 1.4 GHz) 
and from thermal emission. The fraction of thermal emis- 
sion at 1.4 GHz, /therm, was obtained from 
. _ l-(^iM)-^ ... 

^*'-™ - 1 _ {y^lv^)t-l ^^> 

where s is the observed spectral index (in terms of an- 
tenna temperature: Ta oc z/"") between i/i — 1420MHz 
and 1^2 = 408 MHz taken, for each field, from Reich & 
Reich (1988), and 7 = 2.9, t = 2.1 are the spectral in- 
dices of non-thermal and thermal emission, respectively. 
The highest polarization degree of non-thermal emission 
found in these fields is about 20%, well below the theoret- 
ical maximum corresponding to a uniform magnetic field 
(Ginzburg & Syrovatskii 1965): 

"-'^.^ <=) 

which, for 7 = 2.9, is H = 74%. To depolarize to the 
observed level with differential Faraday rotation, RM > 
53 rad m~^ would be required, much higher than the RMs 
directly measured. Thus, other depolarization mechanisms 
must be operating. 



3.2. Distribution of total intensity and polarization 
signals 
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Fig. 3. Distributions of total intensities (isotropic compo- 
nent subtracted, see text) for the U99 regions centered at 
1 = 50° (sofid line), 80° (6 > 8°; short dashes), 146°. 4 (dot- 
ted line), and 200° (long dashes). The maps used for this 
analysis do not include > 5(t unresolved discrete sources; 
addition of these enhances the high intensity tails. 



scope beamwidth is 9.35'), we do not need to worry much 
about that since the derived power spectrum is affected 
only on small scales, where fluctuations due to point 
sources dominate anyway (De Zotti ct al. 1999). 

It is clear from Figs. and that the distributions 
of both total and polarized emissions arc distinctly non- 
Gaussian. The statistical significance of this visual impres- 
sion can be quantified by computing the moments of order 

Mi = ^31 ' ^ ^ (6) 

where Sj is the measured signal (polarized or unpolar- 
ized intensity in the j-th cell) and S is its mean over the 
considered sky region. 

It is convenient to define the dimensionless quantities 
/?! = A^i/A^i: /52 = A*4/M2, Pa = MaA^s/A^i Pi = A^e/A^i, 
/35 = MTA^s/A^i: Pe = m/fJ't- 

Usual definitions of skewness and of kurtosis are 
skewn = V^ and kurt = /?2 — 3, which vanish for a Gaus- 
sian distribution. The probable errors of /3i and /?2 are 
given by (Pearson 1924): 

A/3i = 0.6745(/3i/n)i/2 . 

• (4/34 -24/32 + 36 + 9/3i/32- 12/33 + 35/3i)'/' (7) 



All maps were put onto the same 12' x 12' grid. Al- 
though these cells are not totally independent (the tele- A/32 = 0.6745ri^^" • 
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Table 2. Moments of the distributions of total and polarized intensities. 





mean 


(mK) 


a (mK) 


skewness 


kurtosis 




int 


pol 


int 


pol 


int 


pol 


int 


pol 


Uyaniker 1 = 50° 


2480 


155 


550 


70 


0.9 ±0.5 


1.3 ±0.1 


5.5 ±3.5 


3.2 ±0.3 


no sources 


2360 




440 




0.051 ±0.007 




-0.89 ±0.02 




Uyaniker 1 = 80°* 


1640 


125 


310 


60 


19 ±12 


0.54 ± 0.05 


74 ±36 


1.0 ±0.1 


no sources 


1610 




270 




2.0 ±0.06 




2.0 ±0.1 




Uyaniker 1= 146°. 4 


1840 


110 


190 


65 


28 ±7 


0.43 ± 0.04 


58 ±10 


-0.07 ±0.08 


no sources 


1740 




50 




0.67 ±0.07 




1.0 ±0.1 




Uyaniker 1 = 200° 


1270 


70 


230 


40 


23 ±11 


0.36 ± 0.03 


82 ±22 


-0.05 ±0.1 


no sources 


1200 




140 




0.23 ±0.01 




-0.94 ±0.03 




D97 240° < 1 < 260° 


161 


13 


110 


9 


60 ±40 


0.26 ±0.02 


210 ± 80 


0.08 ±0.07 


no sources 


129 




70 




0.004 ± 0.003 




-0.15 ±0.03 




D97 280° < 1 < 300° 


425 


106 


1741 


33 


340 ± 40 


0.98 ± 0.05 


470 ± 60 


1.0 ±0.1 


no sources 


192 




227 




1.41 ±0.04 




0.53 ±0.07 




D97 300° < 1 < 320° 


577 


15 


930 


15 


90 ±20 


0.88 ±0.05 


200 ± 30 


0.7 ±0.2 


no sources 


423 




404 




0.59 ± 0.02 




-0.22 ±0.04 




D99 20° < 1 < 30° 


317 


32 


907 


20 


130 ± 40 


1.3 ±0.3 


270 ± 70 


4±2 


D99 55° < 1 < 65° 


72 


32 


18 


11 


1500 ± 300 


0.49 ± 0.03 


2400 ± 300 


0.6 ±0.1 



Analysis limited to the region 8° < 6 < 15° 



100 r 




The deviations from the Gaussian value (zero) of skew- 
ness and/or kurtosis are in general highly significant, both 
for total and for polarized emissions. This conclusion is 
strengthened by the Monte Carlo simulations reported in 
Fig. ^ for 10° X 10° patches (grid 4' x 4' as default in the 
D97 data), showing that the skewness differs from at a 
<^ lOcr level in the case of polarization and even more for 
total intensity. This may be a difficulty for methods in- 
volving Wiener filtering of the data to remove foreground 
contributions to CMB maps (Bouchet et al. 1999) since a 
Gaussian approximation for the distribution of foreground 
signals is assumed. On the other hand, Independent Com- 
ponent Analysis algorithms (Baccigalupi et al. 2000) re- 
quire that all independent components contributing to the 
observed maps, except, at most, one have non-Gaussian 
distributions. 



Fig. 4. Global skewness of the total (upper, solid line) and 
polarized emissions (dashed) for the D97 survey compared 
with the rms error on the skewness derived from Monte 
Carlo simulations for a 10° x 10° patch (grid 4' x 4' as 
default in the D97 data) assuming the same power spec- 
trum for the signals but a Gaussian distribution (straight 
dotted line at the bottom). 



•(/?6 - 4/32/?4 + 4/33 - Pi + 16/3i/32 - 8/33 + 16/?i)i/2. (8) 
The values of skewness and kurtosis for the regions 
we have analyzed, and their errors, are given in Table 0, 
together with the mean total and polarization intensities 
and their standard deviations a. Note that, for the U99 
area centered at 1 = 80°, we have limited our analysis to 
the region above 6 = 8° to avoid the relatively large hole 
in the map at lower Galactic latitudes. 



4. Intensity and polarization power spectra 

It is currently standard practice to adopt as a statisti- 
cal measure of the temperature pattern on the celestial 
sphere, the power spectrum of temperature fluctuations, 
Ci, which turned out to be of fundamental importance in 
studies of CMB anisotropics. 

The Cf's are defined as follows. The spherical harmonic 
expansion of the sky signal s writes: 

s{h)=J2^irnYr{n), (9) 

em 
where aim are the expansion coefficients, h is the 

unit vector in the direction (??, 0), and the spherical har- 
monic Yf^ is related to the Legendre function, PJ™ , 
of degree i and order to (|?7i| < £) by Y™{d,(t)) = 

c^Pf^ {cos -d) exp(JTO0); the coefficients c™ normalize the 
spherical harmonics to / dilY™Y™ = SuiSmm! . The cor- 
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Fig. 5. Reconstructed power spectrum of CMB tempera- 
ture fluctuations from a simulated 10° x 10° patch (solid 
line) compared with the theoretical spectrum predicted 
by a standard CDM model, averaged over the whole sky 
(dashed line). 



relation function C{9) =< s{n)s{n') >hfi'=cose can be ex- 
panded into Legendre polynomials, with coefficients given 

by 



Ce 



2£ 



1 " 

— y 

+1 ^^, 



\aen 



(10) 



The multipole £ corresponds to the angular scale 
6* ~ 180/^ degrees . (11) 

We have projected each sky patch onto a null signal 
sphere and computed the power spectrum using of the 
HEALPix tools (Gorski et al. 1998). The coefficients were 
renormalized by simply dividing by the fractional coverage 
of the sky. 

Note that the use of spherical harmonics is not strictly 
necessary, given that we are dealing with limited areas 
of the sky. In fact, for the common data sets, our results 
agree with those by Tucci et al. (2000) who resorted to 
a standard Fourier analysis technique. We preferred, how- 
ever, to stick to the spherical harmonic analysis to be used 
for the all sky MAP and Planck data. 

We have tested our method using a simulation of a 
purely cosmological signal. We have generated an all sky 
map of the CMB temperature distribution as predicted by 
a standard model for a Cold Dark Matter (CDM) dom- 
inated universe, at a resolution of about 3.5'. We have 
then applied our algorithm to a randomly chosen 10° x 10° 
patch. In Fig. ^ we show the recovered power spectrum for 
the patch compared with the theoretical one (which is, of 
course, an all-sky average). The reconstructed spectrum 



oscillates around the theoretical one, due to the sample 
variance, i.e. to the very limited sampling of a random, all 
sky process. However, the main features of the spectrum 
are recovered. 

Linear polarization is described by the Stokes param- 
eters Q and U, from which the total polarization inten- 
sity PI = y/O^+lP and the polarization angle (j)p = 
arctan {U/Q), can be derived. PI is a scalar quantity that 
can be expanded into spherical harmonics [Eq. (^]; its 
power spectrum coefficients Cf can be computed as in 
Eq. i^. 

In order to derive the true power spectrum coefficients 
from the measured ones, C'^'^^ , we must allow for the con- 
tribution of instrumental noise and for the effect of the 
detector response function b{£): 



C, = 



/-^niap 



/^^noisc 



W{£) 



W{£) = b{£y 



(12) 



The beam response for each multipole is modelled by 
6(^)=expM(£-|-l)aLam/2]. (13) 

where CTboam = FWHM/(2V2hr2) and FWHM is the 
full width at half maximum, in radians. In the case of 
D97, the beam is slightly elliptical and Eq. ( p^ ) has been 
modified accordingly. 

The instrumental noise has an approximately flat 
power spectrum: 

2 



C; 



N, 



pixc. 



Is 



(14) 



For the surveys of D97, D99 and U99 we have adopted 
the values of the instrumental noise given by the authors 
(taking into account the different sensitivities reached by 
D97 in different regions). In the case of the BS76 data, we 
found it convenient to treat the noise as a parameter to 
be determined (see subsection 4.2). 

We have analyzed 10° x 10° patches of the surveys by 
D97, D99 and U99. Correspondingly, the minimum value 
of £ for which the power spectrum can be estimated is ~ 
100, corresponding to an angular scale of about 2°: as the 
angular scale approaches the size of the patch, the effects 
of poor sampling become unacceptably large. Information 
on larger angular scales is provided by the BS76 data. 
The maximum value of £ is determined by the angular 
resolution of the survey; we have ^max — 800. 

4-1. Low and medium Galactic latitudes 

We have focused our analysis on the low RM regions of the 
D97, D99, and U99 surveys. As detailed in the following, 
we find that these regions do possess remarkably similar 
polarization power spectra. 

Following Maino et al. (1999) we divided the D97 data 
in twelve 10° x 10° subpatches, and we evaluated the power 
spectrum of both total intensity and polarization for each 
patch, over the range 100 < £ < 800. The results are 
plotted in Figs. ^-g. In each panel, the dashed line shows 
a power law fit {Cg = a ■ £^). 
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Fig. 6. Polarization (upper panels) and total intensity Fig. 8. Polarization (upper panels) and total intensity 

(lower panels) angular power spectra for the D97 region (lower panels) angular power spectra for the D97 region 

360° > 1 > 320°. The panels correspond to 10° intervals in 280° > 1 > 240°. The panels correspond to 10° intervals in 

longitude, with 1 decreasing from left to right. The dashed longitude, with 1 decreasing from left to right. The dashed 

lines represent power law fits. lines represent power law fits. 




Fig. 7. Polarization (upper panels) and total intensity 
(lower panels) angular power spectra for the D97 region 
320° > 1 > 280°. The panels correspond to 10° intervals in 
longitude, with 1 decreasing from left to right. The dashed 
lines represent power law fits. 

Fig. shows the fluctuations around the mean of both 
total and polarized intensities for the D97 survey, as a 
function of the angular distance from the Galactic cen- 
ter, for regions 10° wide in longitude. It may be noted 
that polarization and total intensity fluctuations are cor- 
related only near the Galactic center. The intensity peaks 
are associated to the most prominent Galactic sources, 
i.e. the Galactic center and the Vela Supernova rem- 
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280 300 320 340 360 



Fig. 9. Total (upper panel) and polarized intensity fluc- 
tuations in the area surveyed by D97. 



nant (260° < 1 < 270°). For polarization, the signals at 
320° < 1 < 360° are due to features appearing near the 
Galactic center such as the polarization "plume" (Dun- 
can et al. 1998), while away from the Galactic center, at 
1 < 320° , only the supernova remnant produces peaks both 
in total and in polarized intensity. 
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Fig. 10. Power spectrum for total intensity (dotted line) 
and polarization (dashed line) in the 240° < 1 < 250° 
D97 region. The solid line shows the power spectrum co- 
efficients Ci of total intensity before (dotted line) and 
after (solid line) removing the most prominent sources 
as described in the text, in comparison with polarization 
(dashed line). 



Table 3. Bright Hii regions between 1 = 300° 
320°. 



and 1 = 



1 (deg.) 


h (deg.) 


Sa 


7GHz (Jy) 


Size (') 


305.1 


0.1 




16.3 


5 


305.2 


0.0 




62.2 


7 


305.2 


0.2 




50.1 


4.8 


305.4 


0.2 




62.2 


3.5 


305.6 


0.0 




37.6 


8 


307.6 


-0.3 




12.2 


4.2 


308.6 


0.6 




17.8 


7.8 


308.7 


0.1 




12.0 


- 


308.7 


0.6 




21.0 


- 


309.6 


1.7 




51.0 


7.5 


310.8 


-0.4 




16.6 


10.9 


311.9 


0.1 




12.9 


4 


311.9 


0.2 




11.5 


4 


312.3 


-0.3 




10.0 


- 


316.3 


0.0 




12.0 


- 


316.8 


-0.1 




43.2 


2.7 


317.0 


0.3 




15.0 


7 


318.0 


-0.8 




11.0 


10.8 


319.2 


-0.4 




12.4 


5.9 


320.2 


0.8 




11.0 


1.8 


320.3 


-1.4 




12.0 


- 


320.3 


-1.0 




23.0 


- 


320.4 


-1.1 




17.5 


6 


320.4 


-1.0 




13.0 


5.6 



Also, and most important, this figure shows that the 
total intensity fluctuations are higher by 2 or 3 orders of 
magnitude than the polarization ones, implying a polar- 
ization degree much lower than the maximum expected 
for undepolarized synchrotron emission. 

In fact, most of the observed intensity close to the 
Galactic plane appears to come from bright Hii regions. 
For example, the region between 1 — 300° and 1 = 
320° contains 95 catalogued bright Hii regions (Paladini 
et al. 2000); the brightest ones are listed in Tabled with 
their flux density at 2.7 GHz and their angular sizes which 
are typically in the range 2'-10'. It is easily checked that 
these sources account for most, if not all, of the total in- 
tensity reported by D97. The intensity fluctuations due to 
them can be roughly estimated to be: 

""HII - j:f— > (15) 

^ ^ pixels 

where Si is the flux of the i — th Hii region in the 
catalogue, scaled to the D97 resolution and to 2.4 GHz by 
assuming a spectral index of 0.1 {Si, ex i/""'^), appropriate 
for the free-free emission. This yields ana — 1 Jy/beam, 
close to the value of 1.4 Jy/beam that we obtain from the 
D97 data. 

The free- free emission from Hii regions is not polarized, 
but Thomson scattering by the electrons in the Hii region 
itself may polarize the radiation tangentially to the edges 
of the cloud structure at a maximum level of about 10% 
for an optically thick cloud (Keating et al. 1998). 

Another indication that sources contributing most of 
the intensity in the area surveyed by D97 are unpolarized 
is obtained by removing the highest intensity peaks and 
replacing their signal with the median value in an annulus 
around them. As shown by Fig. |l^, this has the effect of 
strongly decreasing the amplitude of the total intensity 
power spectrum, particularly at intermediate and small 
angular scales, making its shape quite similar to that of the 
polarization power spectrum. Also, the polarization degree 
becomes ~ 0.3, consistent with moderately depolarized 
synchrotron emission. 

As already mentioned, the power spectra obtained for 
each patch have been fitted with simple power laws: 

Ci = a£>^ 100 < £ < 800 . (16) 

The values of the parameters a and /3 have been ob- 
tained by minimizing the quantity: 



'Rt 



Y^ [logio Q - logio ' 



/31ogio^] 



(17) 



i=e„ 



The best-fit values of the parameters are listed in Ta- 
bles H (for total intensity) and ra for polarization. The Ta- 
bles also give the errors Aa, A/3 calculated following the 
standard prescriptions [see, e.g. Press et al. (1996), Chap- 
ter 151, with 



ACf 



E 



(logio Ci - logi 



/31ogio^)' 



L 



(18) 
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Table 4. Power spectrum parameters of total intensity 
for the D97 survey. The first column gives the Galactic 
longitude of the patch center. 



1 (deg.) Q (K^) (3 Aa (K^ 



A/3 



355 


3.1 


-1.66 


0.1 


345 


6.8-10-2 


-1.32 


4-10-^ 


335 


7.3 


-2.04 


0.7 


325 


0.21 


-1.63 


2-10-2 


315 


1.0 


-2.06 


0.2 


305 


0.39 


-1.67 


2-10^2 


295 


0.16 


-1.37 


1-10-2 


285 


11.6 


-1.94 


0.7 


275 


1.M0-* 


-0.76 


5-10-" 


265 


4.2-10-2 


-1.19 


2-10-^ 


255 


7.6-10-^ 


-1.87 


8-10-^ 


245 


5.3-10-5 


-0.98 


4-10"'' 



-10 
-10 
-10 
-10 
-10 
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Table 5. Power spectrum parameters of polarized inten- 
sity for the D97 survey. The first column gives the Galactic 
longitude of the patch center. 



1 (deg.) 


a(K^) 


/3 


Aa (K^) 


A/3 


355 


7.510-" 


-1.49 


7-10-" 


2-10-^ 


345 


2.1 


10-5 


-1.47 


2 


10-** 


2 


10-2 


335 


4.2 


10-5 


-1.44 


4 


10-** 


1 


10-2 


325 


8.1 


10-^ 


-1.93 


8 


10-5 


2 


10-2 


315 


5.2 


10-5 


-1.64 


4 


10-" 


1 


10-2 


305 


3.3 


10-'^ 


-1.28 


4 


10-'^ 


2 


10-2 


295 


3.9 


10-5 


-1.77 


4 


10-" 


2 


10-2 


285 


1.1 


10-^ 


-1.91 


1 


10-5 


2 


10-2 


275 


3.1 


10-5 


-1.67 


4 


10-" 


2 


10-2 


265 


1.5 


10-^ 


-1.43 


3 


10-5 


3 


10-2 


250 


4.0 


10-5 


-1.86 


5 


10-" 


2 


10-2 


245 


3.1 


10-5 


-1.86 


5 


10-" 


2 


10-2 



The visual impression that the fit is generally very 
good is quantitatively confirmed by the fact that in all 
cases Aa <C a and A/3 <C f3. 

In the case of total intensity, the slope /3 shows large 
variations (from ~ —0.8 to ~ —2) along the Galactic 
plane. Steeper slopes correspond to regions where diffuse 
emission dominates; point sources add power on small 
scales (large i). The effect of point sources accounts for 
our finding of total intensity power spectra flatter (and, 
in some regions, much flatter) than derived from previ- 
ous analyses of lower resolution surveys. In fact, analy- 
ses of the 408 MHz Haslam et al. (1982) map (Tegmark 
& Efstathiou 1996) and of the 1420 MHz Reich & Re- 
ich (1988) map (Bouchet & Gispert 1999) yielded values 
of /3 ~ —3. Similar values of /3 are also found from U99 
"background" maps (i.e. those obtained removing > 5a un- 
resolved sources); in this case, we find /3 ~ —2.6, —2.8, 
and —3.35 for the patches centered at 1 = 50°, 1 — 146°. 4 
1 — 196°. 8, respectively. The situation is quite different for 
polarization. In general, the variations of /3 are substan- 



Table 6. Power spectrum parameters for total intensity in 
the 2.7 GHz D99 data. The first column gives the Galactic 
longitude of the the patch center. 



1 (deg.) a (K2) /3 Aa {K'-') Af5 




Table 7. Power spectrum parameters of the polarized in- 
tensity in the 2.7 GHz D99 data. The first column gives 
the Galactic longitude of the patch center. 




tially smaller. Also, there are several regions (at 1 < 320°, 
with exception of the Vela supernova remnant) where the 
power spectra are remarkably similar and described, for 
100 <i< 800, by: 



C[ ^ (l.l±0.7)-10-^-( — 



-1.7±0.2 



V 2.4GHz / 



-5.8 



K^ , (19) 

.2.4GHzy .... . 

where we have explicitly indicated the typical fre- 
quency dependence of Galactic synchrotron emission. We 
find essentially the same power spectra also for several 
D99 regions and for the U99 regions with low rotation 
measures (see below). Therefore, we consider it as typical 
of the diffuse Galactic polarized synchrotron emission. 

The D99 survey at 2.7 GHz covers a region close to 
the Galactic plane at -5° < 6 < 5°, 4.9° < 1 < 74°. We 
selected five squared 10° x 10° patches centered at 1 = 25°, 
35°, 45°, 55°, and 60°. 

The angular power spectra are reported in Figure ^ 
and the best fit values for the coefficients of the power 
law are in Tables o and M. The dotted lines plotted in the 
panels showing the polarization power spectra represent 
Eq. (|l^). Clearly, the polarization power spectra derived 
from the D99 survey are remarkably close to those found 
for the D97 survey, scaled using a typical synchrotron 
spectral index of 2.9. 

Let us consider now the three patches from the U99 
survey at 1.4 GHz having low rotation measures, with cen- 
ters and amplitudes given by: 



1 = 50°, 6=12°, Al= 10°, A6 = 15°, 
1 = 146.4°, b = 6.8°, Al = 12°, Ab = 6.3° 
1=196.8°, 6=11.2°, A1 = A6= 6.7°. 



(20) 
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Fig. 11. Polarized (top) and total intensity (bottom) an- 
gular power spectra for the D99 regions 20° > 1 > 30°, 
30° > 1 > 40°, 40° > 1 > 50°, 50° > 1 > 55°, 55° > 1 > 65° 
from left to right respectively. The dashed lines represent 
power law fits. The dotted hne represents Eq. (18). For 
each panel, I runs from 100 to 1000 on a logarithmic scale. 

Table 8. Power spectrum parameters of polarized (PI) 
and total (I) intensity for the U99 survey at 1.4 GHz. 



region 


a(K^) 


/? 


Aa (K") 


A/3 


1(PI) 


5.6-10-* 


-1.48 


8-10-^ 


2-10~^ 


2 (PI) 


5.3-10-^ 


-2.46 


5-10"^ 
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3 (PI) 


1.9-10"^ 


-2.27 


2-10"^ 


2-10"^ 


1(1) 
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-1.02 


MO"* 


3-10-2 


2(1) 
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-0.51 


7-10-^ 


3-10-2 


3(1) 
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-0.94 
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In Fig. lU we plot the angular power spectra of these 
regions, for polarized (top) and total (bottom) intensity. 
As in the previous figures, the dashed lines show the power 
law fits with the values of the parameters listed in Table ||, 
and the dotted lines represent Eq. ( p^ ) . Again there is a re- 
markable agreement with the power spectra derived from 
the D97 survey, indicating that the polarized synchrotron 
emission keeps essentially the same power spectrum (both 
in amplitude and in slope) up to 6 ~ 10°. 

If we take into account also the power spectra obtained 
for the U99 regions listed above and for the D99 regions 
at 1 > 50°, the average power spectrum is described by: 



C 



(1.2 ±0.8) -10" 



^450 



'1.8±0.3 



V 2.4GHz y 



.2.4GHzy ^^ ' (21) 

very close to that found from the D97 survey only 

[Eq. (©] 



Fig. 12. Polarized (top) and total (bottom) intensity an- 
gular power spectra for the three U99 areas at 1.4 GHz 
with low rotation measures specified by Eq. (19): regions 
1, 2, 3 from left to right. The dashed lines represent power 
law fits. The dotted line shows the average power spec- 
trum for the D97 survey, scaled in frequency according to 
Eq. (18). For each panel, i runs from 100 to 1000 on a 
logarithmic scale. 



^.2. High Galactic latitudes 

The accurate and high resolution D97, D99 and U99 polar- 
ization surveys allow to understand the correlation proper- 
ties of polarized synchrotron emission at multipoles larger 
than ~ 100, owing to the limited extent of the patches we 
can extract from the observed regions, and cover regions 
at low or intermediate Galactic latitudes. 

In order to extend the analysis of the polarized syn- 
chrotron fluctuation power spectrum to superdegree angu- 
lar scales and to high Galactic latitudes we have exploited 
the BS76 polarization measurements. 

We restrict the present analysis to the 1.411 GHz chan- 
nel, for a simpler comparison with the high resolution 
surveys. As we already mentioned in the Introduction, at 
1.411 GHz the BS76 measurements refer to 1726 positions 
over approximately half of the sky at 1 /cL 0° with a beam 
FWHM of about 0.6°. However, the sky is significantly 
undersampled and the coverage is inhomogeneous. 

We projected the original data into HEALPix maps 
at different resolutions 9 ~ 3600/nsido arcminutes, 
where Uside is the HEALPix resolution parameter (Gorski 
et al. 1998). The whole map is essentially filled only at low 
resolutions, 9 ~ 7.3° or 3.6°, corresponding to Uside = 8 or 
16 respectively, that allow to estimate the angular power 
spectrum only at ^ ;S 50. On the other hand, we can take 

advantage of some regions where the sky is better sampled 
and/or interpolate the existing data to fill maps at higher 
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resolutions, 9 ~ 1.8°-0.9°, corresponding to riside = 32 
or 64 respectively, to try an approximate estimate of the 
power spectrum at higher multipoles, up to ^ ~ 100-200, 
and reach the multipole range where the power spectrum 
estimation from high resolution surveys does not suffer of 
significant boundary effects introduced by patch sizes. It 
is important in any case to fill the map through interpo- 
lations to avoid the "holes" corresponding to unobserved 
pixels. Such holes would introduce spurious (flat) power 
on the pixel scale, being seen like "negative sources" in 
the power spectrum estimation (see La Porta & Burig- 
ana 2000 for further details). We do not expect crucial 
artifacts from this treatment, since the analysis of the po- 
larized synchrotron fluctuations on smaller angular scales 
presented above indicates that the power significantly de- 
creases toward high multipoles. Of course, by comparing 
the power spectrum derived from the whole map and from 
the regions where the sky is better sampled (and, conse- 
quently, the possible interpolation effects less relevant) we 
can test the effect of these approximations. 

We identified three relatively well sampled regions in 
the maps derived from the BS76 data: 1) a patch at low 
Galactic latitude (110° < 1 < 160°, 0° < 5 < 20°), that 
can be used for a comparison with the above results at 
higher resolutions and to extend the analysis of the Galac- 
tic plane regions to low multipoles; 2) a relatively wide 
patch around the North Galactic Pole [(5° < 1 < 80°, 
b > 50°) together with (0° < 1 < 5°, 5 > 60°) and 
(335° < 1 < 360°, b > 60°) ]; 3) a relatively smaU but bet- 
ter sampled region (10° < 1 < 80°, b > 70°), included in 
the previous patch. The analysis of the patches 2 and 3 al- 
lows to extend the polarized synchrotron power spectrum 
estimation to high Galactic latitudes, where the Galactic 
contamination is expected to be lower and, correspond- 
ingly, the view of the CMB is cleaner (however, the min- 
imum of the Galactic emission is near 1 = 190°, b = 50°, 
see Berkhuijsen 1971). 

The results, scaled to the case of a full sky coverage, 
and the fits discussed below, are shown in Fig. |l^, for 
SO ^ i ^ 200. At multipoles ^ 30 boundary effects are 
relevant, particularly for the patches. At low £'s the power 
spectrum decreases rather steeply with increasing multi- 
pole number, then, at ^ <L 100, it flattens out. Given the 
sensitivities quoted by BS76, it is likely that the flattening 
is mostly due to instrumental noise. 

To investigate the effect of the latter we have fitted the 
derived power spectra as the sum of a power law plus white 
noise due to the instrument. The best fit values of the 
parameters were computed using the MINUITS software 
package of the CERN libraries (available at the WEB site 
tittp : //cern . web . cern . ch/CERN7| ) . We neglected the ef- 
fect of beam smoothing, which becomes important at mul- 
tipoles higher than those at which the noise power starts 
to be dominant. 



In order to test the stability of the results, we have 
carried out the calculations using two different routines: 
MIGRAD and SIMPLEX. The results are very close to 
each other except for patch 1, as discussed below. For the 
other three cases (whole map and patches 2 and 3) we 
present only the fits produced by MIGRAD, while both 
fits are given for patch 1 (see also Fig. 03). 

MIGRAD yields: 
C[ ~ (9.4 • 10"^ • e-^-^ + 3.2 • 10"^) K^ (22) 

for the whole map. 



4.8-10 



-7\ 



K' 



(23) 
(24) 



c[ ~ (1.5 • r^-i 

for patch 2, 
C[ ~ (1.1 • e-'^-^ + 8.7 • lO"'^) K^ 

for the patch 3, and 
C[ ~ (8.8 • 10"^ • r^-®^ + 1.1 • 10"^) K^ (25) 

for the patch 1. For the latter patch the SIMPLEX fit 
is 
C[ ~ (2.0 • 10^2 • e-^-^ + 2.2 • 10"^) K^ . (26) 

As illustrated by Fig. |l^, the power spectrum of po- 
larized emission for patch 1 shows substantial oscillations, 
which are responsible for the instability of the power law 
fit derived with the two routines. We believe the MIGRAD 
result to be more reliable for at least two reasons: first, it is 
consistent with that derived from the data of U99 for the 
area common to both surveys, for the overlapping range of 
scales {£ -- 200), while the SIMPLEX fit yields a signal a 
factor of several too high; second, the estimated noise level 
implied by the SIMPLEX fit seems to be too low, being 
substantially lower than found for the other two patches 
and even smaller than found for the whole map (in all 
these cases SIMPLEX and MIGRAD results are in close 
agreement), whereas the estimated noise level implied by 
the MIGRAD fit is in agreement also with that derived 
by inverting a patch of noise simulated on the basis of the 
sensitivities quoted by BS76 in the sky directions of the 
considered region. 

If we can ignore the SIMPLEX result for patch 1, as 
due to a numerical instability in the presence of large oscil- 
lations of the power spectrum, we may conclude that the 
power spectra for the three regions agree with each other 
to within a factor of about 2. Remarkably, no significant 
differences are found between the region on the Galactic 
plane and the North polar region. The power law terms in 
Eqs.(p3|)-(|25|) are all consistent with 



C 



(4 ±3) -10"^ 



-2.9±0.2 



V2.4GHzy 



-5.8 



.2.4GHzy ^'' (^^) 

where, to ease the comparison with Eq. (El[), the ampli- 
tude was scaled to 2.4 GHz using the typical synchrotron 
spectrum. 

The power spectrum found for the whole map exhibits 
a slope quite close to that of the three patches considered 
and an amplitude about 3-6 times smaller, mostly due to 
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Fig. 13. Fits of the angular power spectrum of polarized synchrotron emission at 1.411 GHz derived by projecting 
into a map the Brouw & Spoelstra (1976) data for different sky regions. Top left panel: whole map; top right panel: 
patch 1; bottom left panel: patch 2; bottom right panel: patch 3. Note the large oscillations found for the patch 1. The 
smooth solid lines show the MIGRAD fits. The smooth upper solid lines include also the noise contribution, whereas 
the lower ones show the power law approximation of the Galactic synchrotron component only, as derived from the 
fit. SIMPLEX fits are very close to these for patches 2 and 3 and are not shown; on the other hand, we show (dashed 
line in the top right panel) the discrepant SIMPLEX fit for patch 1. See the text for further details. 



the lower polarization degree in the regions far from the 
patches considered here. 

In spite of their poor sky sampling and of all the ap- 
proximations introduced to treat them, the BS76 data in- 
dicate an amplitude of the polarization power spectrum 
in the region of overlap {£ ~ 100) quite close to that 
derived from more recent surveys, although the slope is 
steeper. This suggests that, at least on superdegree angu- 
lar scales, the fluctuations of the diffuse Galactic polarized 
synchrotron emission are almost independent of Galactic 
latitude. 



5. Discussion and conclusions 

We have analyzed the available surveys of diffuse Galac- 
tic polarized emission at GHz frequencies. The most re- 



cent ones cover regions at low and medium Galactic lat- 
itudes with angular resolution ~ 10' or better (Duncan 
et al. 1997, 1999, Uyaniker et al. 1999). Observations at 
high Galactic latitudes have a resolution of 0°.6 (Brouw 
& Spoelstra 1976). 

A particularly bewildering issue is Faraday depolar- 
ization, which may be substantial at the frequencies of 
the data analyzed here (a few GHz) but will be negligi- 
ble at 100 GHz. Space varying Faraday depolarization ef- 
fects may substantially affect the power spectrum of polar- 
ized synchrotron emission in several ways: the amplitude 
is decreased and structure may be created on a variety of 
scales. If so, the validity of extrapolations of the present re- 
sults to map's and Planck's frequencies may be largely 
spoiled. However, our careful analysis of depolarization ef- 
fects has led to the identification of regions where rotation 
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Fig. 14. Comparison of the power spectrum of the CMB 
polarized component predicted by a standard CDM model 
(dashed Hne) with the power spectrum of the polarized 
Galactic synchrotron emission at 100 GHz, as yielded by 
Eqs. (21) and (27) for high and low values of ^, respec- 
tively (solid lines). Also shown are simple minded esti- 
mates made assuming that the synchrotron emission has 
everywhere the maximum theoretical polarization degree 
(~ 75%): the long-dashed and the dotted hues correspond 
to the average power spectrum of the total power syn- 
chrotron emission at high Galactic latitudes (|6| > 30°) 
derived by one of us (C. Burigana; see also Tegmark & 
Efstathiou 1996) from the map of Haslam et al. (1982), 
and to the median power spectrum estimated by Bouchet 
& Gispert (1999; see also Bouchet et al. 1999), respec- 
tively. 



measures of pulsars and extragalactic sources, the polar- 
ization degree and, in some cases, data on the distribution 
of polarization vectors and on the Galactic magnetic field, 
consistently indicate that Faraday depolarization must be 
small. 

The low polarization degree of radio emission close to 
the Galactic plane is interpreted as due to large contribu- 
tions to the observed intensity from unpolarized sources, 
primarily strong Hii regions. Since such sources are con- 
centrated on the Galactic plane, estimates of the power 
spectrum of total intensity at low Galactic latitudes are 
not representative of the spatial distribution of Galactic 
emission far from the plane. 

Since we need to extrapolate the results by a large fac- 
tor in frequency (from a few GHz up to 100 GHz) a sec- 
ond important issue is the spectral index of synchrotron 
emission to be adopted. The average spectral index of the 
antenna temperature between 408 MHz and 7.5 GHz is 



about 2.8 (Platania et al. 1998); it is expected to steepen 
to ~ 3 above ~ 10 GHz as a consequence of the steepen- 
ing of the energy spectrum of cosmic rays above ~ 15 GeV, 
due to energy losses (Banday & Wolfendale 1990, 1991). 
Indications of a steepening were indeed found by Platania 
et al. (1998) already at 7.5 GHz and from their compar- 
ison of the 408 MHz map by Haslam et al. (1982) with 
a preliminary map at 19 GHz (Cottingham 1987; Boughn 
et al. 1990). A further steepening is expected at still higher 
frequencies. Thus, our choice of an average spectral index 
of 2.9 up to 100 GHz is rather on the conservative side: 
it probably leads to an overestimate of the amplitude of 
the synchrotron power spectrum at this frequency. On the 
other hand, it must be kept in mind that there are evi- 
dences (Reich & Reich 1988; Platania et al. 1998) of flatter 
spectral indices at higher Galactic latitudes, where CMB 
maps are expected to be the cleanest. Also, synchrotron 
spectral indices show substantial spatial variations, which 
will yield a high frequency power spectrum significantly 
different from the one estimated from the data analyzed 
here. 

Yet another issue is the extrapolation of our results 
to different regions of the Galaxy, particularly to higher 
Galactic latitudes. In general, we expect that the power 
spectrum of total and polarized emission varies across the 
sky. For example, it is reasonable to expect less small scale 
structure in the general anticenter region because emission 
cells have generally larger angular sizes, being, on average, 
relatively closer. 

It is likely that the mean amplitude of the power spec- 
trum of Galactic emission decreases with increasing Galac- 
tic latitude as a consequence of the decreased emission. 
Also, as pointed out by Davies & Wilkinson (1999), the 
magnetic field pattern may be more ordered at high Galac- 
tic latitudes, resulting in less small scale structure. More- 
over, narrow depolarizing structures, that may be present 
in regions where depolarization is generally small, may 
again increase the amplitude of the the polarization fluc- 
tuations on small scales observed at relatively low frequen- 
cies. Since polarization surveys cover only a limited frac- 
tion of the sky and, moreover, observations at high Galac- 
tic latitudes allow to estimate the polarization power spec- 
trum only on super-degree scales, all these issues remain, 
to a large extent, open. 

However, our results do suggest that the dependence 
on Galactic latitude of the power spectrum of polarized 
synchrotron emission is weak. Polarization fluctuations 
are found to remain rather low even at relatively low 
Galactic latitudes. This is very encouraging since the high- 
est sensitivity polarization maps that will be produced 
by the Planck satellite will cover regions around the 
Ecliptic poles, which are at moderate Galactic latitudes 
(|6|c.27°). 

In spite of all these uncertainties, our analysis sug- 
gests that the Galactic polarization fluctuations are un- 
likely to hinder MAP's and Planck's measurements of 



C. Baccigalupi et al.: Power spectrum of the polarized diffuse Galactic radio emission 



15 



CMB polarization for ^ > 50. This is illustrated by Fig. |4| 
where the expected level of CMB polarization fluctuations 
predicted by a standard CDM model is compared with 
our estimates of the Galactic contamination at 100 GHz 
[Eqs. (^ and (|2^)]. The same figure also shows the aver- 
age power spectrum of the total power synchrotron emis- 
sion at high Galactic latitudes (|&| > 30°) derived by one 
of us (C. Burigana; see also Tegmark & Efstathiou 1996) 
from the map of Haslam et al. (1982), scaled to the maxi- 
mum theoretical polarization degree {^ 75%), correspond- 
ing to a uniform magnetic field (long dashed line). This 
polarization degree is a conservative upper limit, since tur- 
bulent components of the magnetic field decrease the po- 
larization degree. Also, there are evidences of significantly 
lower synchotron emissions over relatively large regions of 
the sky. In Fig. |lj this is illustrated by the dotted line, 
showing the median synchotron power spectrum estimated 
by Bouchet & Gispert (1999); see also Bouchet et al. 1999) 
from the 1420 MHz survey of Reich & Reich (1988), ex- 
trapolated to higher frequencies with a spectral index of 
2.9; again a polarization degree of 75% has been adopted. 
The CMB polarized component is usually decomposed in 
suitable eigenfunctions keeping memory of different kinds 
of cosmological perturbations, generally known as E and 
B modes (see Hu et al. 1998 and references therein). We 
find that the Galactic synchrotron emission discussed here 
contributes almost equally to the two modes, as expected 
since Galactic and CMB signals are completely uncorre- 
lated. 

Finally, the distributions of both total and polarized 
Galactic emissions was shown to be non-Gaussian at a 
high significance level. This may be a problem for meth- 
ods of component separation using Wiener filtering, which 
assumes Gaussian distributions. On the other hand, it al- 
lows the application of Independent Component Analysis 
techniques (Baccigalupi et al. 2000) which just rely on the 
assumption that all but at most one of the components to 
be separated have a non-Gaussian distribution. 
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